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ABSTRACT 

fy| We present a careful and detailed light curve analysis of RR Lyrae stars in the Small 

Magellanic Cloud (SMC) discovered by the Optical Gravitational Lensing Experiment 
(OGLE) project. Out of 536 single mode RR Lyrae stars selected from the database, 
we have investigated the physical properties of 335 'normal looking' RRab stars and 
17 RRc stars that have good quality photometric light curves. We have also been able 
Q ■ to estimate the distance modulus of the cloud which is in good agreement with those 

l—j ' determined from other independent methods. The Fourier decomposition method has 

c/2 , been used to study the basic properties of these variables. Accurate Fourier decom- 

position parameters of 536 RR Lyrae stars in the OGLE-II database are computed. 
Empirical relations between the Fourier parameters and some physical parameters of 
these variables have been used to estimate the physical parameters for the stars from 
J> 1 the Fourier analysis. Further, the Fourier decomposition of the light curves of the SMC 

[ RR Lyrae stars yields their mean physical parameters as: [Fe/H] = -1.56 ±0.25, M = 

O 0.55 ±O.O1M , T Gff = 6404 ±12 K, logL = 1.60 ±0.01 L and M v = 0.78 ±0.02 for 

OO ■ 335 RRab variables and [Fe/H] = -1.90 ± 0.13, M = 0.82 ±0.18 M , T off = 7177 ± 16 

; K, log L = 1.62 ± 0.02 L and M v = 0.76 ± 0.05 for 17 RRc stars. Using the ab- 

f^*) ■ solute magnitude together with the mean magnitude, intensity-weighted mean mag- 

nitude and the phase-weighted mean magnitude of the RR Lyrae stars, the mean 
0^ , distance modulus to the SMC is estimated to be 18.86 ±0.01 mag, 18.83 ±0.01 mag 

and 18.84 ±0.01 mag respectively from the RRab stars. From the RRc stars, the corre- 
sponding distance modulus values are found to be 18.92 ±0.04 mag, 18.89 ±0.04 mag 
and 18.89 ±0.04 mag respectively. Since Fourier analysis is a very powerful tool for the 
' study of the physical properties of the RR Lyrae stars, we emphasize the importance 

^ \ of exploring the reliability of the calculation of Fourier parameters together with the 

uncertainty estimates keeping in view the large collections of photometric light curves 
that will become available from variable star projects of the future. 
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1 INTRODUCTION 

RR Lyrae (RRL) stars have played an important role in de- 
termining the cosmic distance scale in modern astronomy. 
While several studies have used classical Cepheids as pri- 
mary distance indicators, RRL stars have received substan- 
tial attention in solving the distance scale problem. Apart 
from distance determinations, RRL stars are of particular 
importance as a test bed for the theories of stellar and galac- 
tic structure and evolution. They are the tracers of old stellar 
populations of the bulge, disk and halo components that are 

* E-mail: sdeb@physics.du.ac. in 



present everywhere. RRL stars are radially pulsating A-F 
variable stars with periods in the 0.2 - 1.2 day range, and 
amplitude of variation ^ 2 mag. They can be easily identi- 
fied, and play a key role as the cornerstone of the Population 
II distance scale. They are extensively used to determine dis- 
tances to old and sufficiently metal-poor systems, where they 
are commonly found in large numbers. In particular, RRL 
stars are present in globular clusters (GCs) and the dwarf 
galaxies in the neighborhood of the Milky Way (Greco et al. 
2007), and have also been identified in the M31 field (Brown 
et al. 2004, Dolphin et al. 2004), in some M31 companions 
(Pritzl et al. 2005), and in at least four M31 GCs (Clemen- 
tini et al. 2001). Distances to the Large Magellanic Cloud 
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(LMC) for the population II objects are based on the lumi- 
nosity of the RRL stars (Clementini et al. 2002). 

There have been a very few studies to estimate the dis- 
tance scale of the SMC using the RRL stars. Using four RRL 
light curves of the SMC cluster NGC 121, Walker & Mack 
(1988) have estimated the distance modulus of the cluster to 
be 18.86±0.07 mag. On the other hand, using 22 RR Lyrae 
stars surveyed around 1.3 square degrees near the northeast 
arm of the SMC field NGC 361, Smith et al. (1992) obtained 
the distance modulus of the SMC as 18.90±0.16 mag. Also, 
there are other studies of estimating the distance scale of 
the SMC based on the binary star light curves and double 
mode Cepheids. Harries et al. (2003) reported that the dis- 
tance modulus to the SMC is of the order of 18.89±0.14 
mag taking 10 eclipsing binaries in the SMC. By selecting 
40 eclipsing binaries of spectral type O and B in the SMC, 
Hilditch et al. (2005) have derived the fundamental parame- 
ters of the binaries and refined the distance modulus to the 
SMC to 18.91±0.1 mag. Also Kovacs et al. (2000) found the 
distance modulus of the SMC to be 19.05±0.017 mag based 
on the photometric data of double-mode Cepheids from the 
OGLE project. 

The stellar atmospheric parameters of effective temper- 
ature (T e ff) and surface gravity (log g) are of fundamental 
astrophysical importance. They are the prerequisites to any 
detailed abundance analysis and define the physical condi- 
tions in the stellar atmosphere and hence are directly related 
to the physical properties mass (M), radius (R) and luminos- 
ity (L) of the star. In this paper, we present an independent 
analysis of 536 SMC RRL stars discovered by the OGLE 
project (Soszynsky et al. 2002). The OGLE database is a 
very wealthy resource for studying the characteristics of vari- 
able stars in the Galaxy, LMC and SMC. For the first time, 
we make use of the OGLE SMC data to estimate the dis- 
tance scale of the SMC using a large number of well-sampled 
RRL light curves. The road map of the present investigation 
is to perform a Fourier analysis of the RRL stars in order 
to estimate their physical parameters and hence the SMC 
distance scale. 

We employ the Fourier decomposition technique which 
is used extensively to characterize the observed photometric 
light curves of RRL and other types of variables. The accu- 
rate determination of the Fourier coefficients is, therefore, 
an important task. We have performed an independent au- 
tomated Fourier analysis of all the RRL light curves selected 
in this paper by a computer code developed by us. In section 
2 we give a brief description of the database that we use and 
the procedure of removing the outliers from the light curve 
data. We present Fourier decomposition of the light curves 
in section 3. We also describe the use of the unit-lag auto- 
correlation function for finding out the optimal order of the 
fit to the RRL light curves. Section 4 describes the error 
analysis of the Fourier decomposition parameters <f>n and 
Rii. Section 5 describes the calibration of the I band data 
to the V band. In section 6, we describe the various physical 
parameters of the RRLs obtained by using empirical rela- 
tions from the literature. Section 7 describes the distance 
determination of the SMC. Lastly, in section 8, we present 
the conclusions of our study. 
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Figure 1. Examples of light curve data from the OGLE-II 
database (a) raw RRab light curve, (c) raw RRc light curve. In 
panels (b) and (d) the corresponding outlier removed smoothed 
light curves are shown. 

2 THE OBSERVATIONAL DATABASE 

The RRL stars analyzed in the present work were selected 
from the high quality photometric catalog of RRL stars in 
the SMC from OGLE-II database (Soszyfiski et al. 2002, 
hereafter SZ02). The catalog contains B, V and I light curves 
of 58 RRc stars and 478 RRab stars located in the 11 areas 
close to the bar of the SMC in Johnson- Cousins photometric 
system. The target stars selected for the present analysis 
were passed through a multi-pass nonlinear fitting algorithm 
in IDL (Interactive Data Language) which is very efficient 
in removing the outliers from the data set. Points lying 2a 
away from the fit are rejected from the data set so that 
the objects with well sampled and accurate light curves are 
available for the analysis. All the selected targets have evenly 
covered I band light curves with about 100 - 400 data points 
in I band with an internal accuracy of 0.03 - 0.13 mag in I 
band. The B and V light curves in the SZ02 catalog are of 
lower quality and have a limited number of data points (20- 
40). Therefore, we use only the I band data for estimation 
of the Fourier parameters. In Fig. 1(a) and (c) we show the 
typical phased light curve of a RRab and a RRc star in the 
database. The corresponding light curves with outliers ^ 2a 
removed are shown in panels (b) and (d) respectively. 



3 FOURIER DECOMPOSITION OF THE 
LIGHT CURVES 

Fourier decomposition technique is a powerful and robust 
technique for describing the shape of the photometric light 
curves of the RRL and other variable stars. The method 
has also been utilized for variable star classification (Deb & 
Singh 2009), finding the physical parameters of stars like ab- 
solute magnitude, metallicity, effective temperature, gravity, 
helium abundance etc. of RRL and other stars (Kovacs & 
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Figure 2. Histogram plot of unit-lag auto-correlation function 
for 536 SMC RRL stars in the OGLE database. The order of the 
fit is 4 and 5 respectively for the RRc and RRab stars. 
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Figure 3. Histogram plot of \v 2 of 536 RRL stars selected from 
the database. 



Kupi 2007 and references therein). The method in its mod- 
ern form has been revived by Simon & Lee (1981) to describe 
the progression of Cepheid light curve with increasing period 
called Hertzsprung progression. They showed that the lower 
order Fourier parameters are able to completely describe 
the progression of Cepheid light curves. With theoretical 
light curves based on hydrodynamical models Simon & Lee 
(1981), Simon & Clement (1993, hereafter SC93), Clement 
& Rowe (2000) were able to derive various physical parame- 
ters of stars. On the other hand Kovacs and his collaborators 
(Jurcsik & Kovacs 1996, hereafter JK96, Kovacs & Jurcsik 
1996, Jurcsik 1998, Kovacs 1998 , Kovacs & Kanbur 1998, 
Kovacs & Walker 1999) were able to estimate astrophysical 
parameters by establishing empirical relationships between 
the Fourier parameters and the physical parameters of stars 
in the field. A number of studies have made use of these 
empirical relations to derive the physical parameters of the 
RRL stars. To list a few of these studies, Kaluzny et al. 
(2000) used these empirical relations to estimate physical 
parameters of 26 RRab and 16 RRc stars in M5, Pena et al. 
(2007) used these empirical relations to estimate physical pa- 
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Figure 4. Fourier fitted light curves of a selection of RRab stars 
after pre-processing of outliers removal. The order of the fit to 
the light curve is 5. The solid lines show the Fourier fitted light 
curves. 



rameters of 7 RRL stars in Bootes and Arellano et al. (2008) 
have determined parameters for few RRL variables in NGC 
5466. There are a number of other studies that have made 
use of the Fourier coefficients to determine various param- 
eters of astrophysical importance. Sandage (2004) studied 
the correlation of particular Fourier components of the light 
curves of RRab stars with metallicty discovered by Simon 
and later by Kovacs and his co-workers. On the other hand, 
Morgan et al. (2006, 2007) have derived empirical relations 
for RRc stars connecting [Fe/H], period (P) and Fourier pa- 
rameter 03i. In all the above studies, it has been seen that 
Fourier parameters can indeed be linked with global stellar 
parameters such as luminosity, mass, temperature, metal- 
licty and radius. Most of the physical parameter estimations 
of stars require accurate calculation of Fourier parameters 
and any inaccuracies in the determination of Fourier param- 
eters will reflect in the calculated physical parameters of the 
stars. 

We computed Fourier decomposition coefficients for 536 
RRL stars in the SMC for the I band photometric data. The 
observed magnitudes were fitted with a Fourier cosine series 
of the form 



I(t) = A + V] A, cos(iiu(t - t ) + 4>i), 



(1) 
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Figure 5. Fourier fitted light curves of a selection of RRc stars 
after removing the outliers. The order of the fit to the light curve 
is 4. The solid lines are the Fourier fitted light curves. 



where I(t) is the observed magnitude, Ao is the mean mag- 
nitude, uj=2n/P is the angular frequency, P is the period of 
the star in days, t is the time of observation, to is the epoch 
of maximum light, Ai and <f>i are the zth order Fourier co- 
efficients and m is the order of the fit. Eqn. (1) has 2m + 1 
unknown parameters which require at least the same num- 
ber of data points to solve for these parameters. 

Since period is known from the database, the observa- 
tion time can be folded into phase ($) as 



9= <t-toL_ Int ({t-t ) 



Here to is the epoch of maximum light of the RRL light 
curves. The value of <& is from to 1, corresponding to a full 
cycle of pulsation and and Int denotes the integer part of the 
quantity. Hence, Eqn. (1) can be written as (Schaltenbrand 
& Tammann 1971) 



I(t) = A + A t cos{2Tri$(t) + * 



(2) 
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Figure 6. Fourier phase differences cfel and 03i as a function 
of log P for the I band data of 536 program stars. Open circles 
represent RRab stars and filled circles RRc stars. 



tween and 27t so that they are comparable to those avail- 
able in the literature. 

The optimal order of the fit (m) was chosen to be 4 for 
RRc stars and 5 for RRab stars by the calculation of unit- 
lag auto-correlation function (Fig. 2) and looking at the \ 2 
value of the fit. A detailed discussion on the optimal order of 
the fit has been given in Deb & Singh (2009) which involves 
the calculation of unit-lag auto-correlation function using 
Baart's condition (Baart 1982, Petersen 1986). Increasing 
the order of the fit may reduce the % 2 to some extent but this 
will underestimate the distribution. The fits to the phase- 
folded light curves were made using Levenberg-Marquardt 
(LM) algorithm which is based on x 2 minimization method 
(Press et al. 1992). The reduced chi-square for the fit to the 
RRL light curves are nearly 1. The histogram plot of the x 2 „ 
of the RRL stars is shown in Fig. 3, where v is the degree 
of freedom of the fit which is defined as the number of data 
points minus the number of parameters used to fit the light 
curve data. 

In Figs. 4 & 5 we show smoothed (outliers removed) 
Fourier fitted light curves for 16 randomly selected RRab 
and RRc stars respectively. Further, the Fourier parameters 
for all the stars computed from the I band data are given in 
Table 1 (for RRab) and Table 2 (for RRc). N obs denotes the 
number of data points for each of the light curves retained for 
the analysis after removing the outliers. D m is the deviation 
parameter for RRab stars and is discussed in Sec. (6.1). In 
Fig. 6 we plot the Fourier amplitude ratios R21, R31, R41 
and phase differences ^21 , <fei , 041 versus log P for the RRL 
dataset. 



The Fourier parameters are defined as 

Rii = -p; 4>n = 4>i — i<f>i, 
Ai 

where i > 1. The (/>ji values have been adjusted to lie be- 



4 ERRORS IN THE FOURIER PARAMETERS 

The error in x = f(u, v), where the standard errors in u's 
and v's are known from LM method, is given by 
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Table 1. Fourier parameters for 478 fundamental mode (RRab) SMC RRL variables in the OGLE database (I BAND DATA) 
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Complete table is available in the electronic form. 



Table 2. Fourier parameters for 58 overtone mode (RRc) SMC RRL variables in the OGLE database (I BAND DATA) 
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0.016 


0.123 


0.136 


0.135 


1.416 


4.446 


3.286 


004631.02-724658. 


8 


0.308220 


300 


0.459 


0.059 


0.362 


19.241 


0.011 


0.322 


0.096 


0.071 


4.444 


2.469 


2.491 


003946.61-733332. 





0.309804 


273 


0.514 


0.062 


0.200 


19.420 


0.008 


0.046 


0.099 


0.080 


4.162 


1.948 


3.494 


004902.13-724513. 





0.310884 


314 


0.599 


0.057 


0.282 


19.180 


0.005 


0.541 


0.157 


0.053 


4.196 


2.315 


6.005 


004245.08-730939. 


7 


0.311416 


286 


0.739 


0.060 


0.339 


19.277 


0.016 


0.297 


0.031 


0.104 


4.817 


3.105 


1.882 


005503.23-731907. 


8 


0.313536 


278 


1.096 


0.061 


0.321 


19.275 


0.011 


0.171 


0.076 


0.075 


4.616 


2.970 


5.392 


005913.99-724719. 


4 


0.313798 


252 


0.584 


0.064 


0.305 


19.440 


0.019 


0.078 


0.083 


0.128 


6.105 


3.581 


5.656 


010249.92-720947. 


1 


0.316409 


277 


0.518 


0.061 


0.271 


19.379 


0.014 


0.173 


0.039 


0.110 


4.034 


1.202 


2.223 


004330.74-733416. 





0.320723 


2X11 


0.741 


0.060 


0.174 


19.293 


0.002 


0.096 


0.115 


0.021 


5.961 


2.734 


1.187 



Complete table is available in the electronic form. 



2, s 2,dx 2 ,dx 2 2 dx dx 

a ( x)=(Tu (_) +av (_) + + 2a uv (^)(^)-(3) 

In the case of a large number of observations Eqn. (3) can 
be reasonably approximated by 



2/ \ 2,dx. 2 2(9x^2 

a (x) = a u — ) +a v — ) . 

ou av 



(4) 



Therefore the errors in and <f>n can be approximately 
written as 



^yVi^ + ^Ai 2 ), 



(5) 



o>i + i o>i' 



(6) 



In Figs. 7 & 8 we plot the distribution of estimated errors in 
the quantities R21, R31, R41 and <f>2i, 4>3i and 4>4,i- As is quite 
obvious that the errors in <f>si are larger than the errors in 
021 because the amplitudes for the higher order coefficients 
are smaller and hence it is very difficult to derive their phase 
with better precision. 
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Figure 7. Histogram plot of the distribution of standard errors 
in R21 , R31 and R41 for the I band data. 
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Figure 8. Histogram plot of the distribution of standard errors 
in 021, </>3i and 04ifor the I band data. 



5 CALIBRATION OF THE OGLE SMC RRL 
DATASETS 

5.1 Fourier parameters in V band 

The Fourier coefficients of a particular light curve in different 
photometric bands are different. Nonlinear hydrodynamical 
models have been used by Dorfi & Feuchtinger (1999, here- 
after DF99) to obtain UBVI light curves of RRL stars. They 
have computed the inter-relations of the Fourier parameters 
in V and I bands and compared them with those obtained 
by Morgan et al. (1998, hereafter MSB98) from observed 
light curves of metal poor globular cluster M68 by Walker 





1 9.0 
mag (A ) 



Figure 9. Histogram plot of magnitudes of 536 RRL stars. The 
solid line is the best fit normal distribution. 



(1994, hereafter W94). However, only 8 RRab and 8 RRc 
stars from W94 were used for setting up the relations in 
MSB98. While the V-I inter-relations for amplitude ratios 
(Rn) of DF99 showed good agreement with the empirical 
relations of MSB98, large deviations showed up in the rela- 
tions of some of the phase differences (4>n) for two sets of 
models with Z = 0.001 and 0.0001. DF99 concluded that no 
distinct metallicity effect could be followed from these re- 
sults and emphasized the need for more model calculations 
covering a wide range of stellar parameters. 

In order to check the reliability of the theoretical inter- 
relations of the phase parameters of DF99, we have calcu- 
lated the Fourier phase parameters in the I and V bands 
from the highly accurate photometric light curves of selected 
RRab and RRc stars in globular cluster M3 from Benko et al. 
(2006, hereafter B06) together with the data of 8 RRab and 
8 RRc stars from W94. All the RRab stars free from Blazhko 
effect and RRc stars having good quality light curves were 
selected from B06 after visual inspection. We have used 29 
RRab variables 1, 6, 9, 25, 27, 31, 32, 42, 53, 57, 58, 69, 84, 
89, 93, 94, 100,109,135,137, 139, 142, 144, 146, 148, 165, 167, 
175, 222 of globular cluster M3 from B06 and 8 variables 2, 
9, 10, 14, 22, 23, 25, 35 of globular cluster M68 from W94. In 
the case of RRc stars, 16 variables 21, 37, 75, 86, 88, 107, 126, 
128, 131, 147, 152, 171, 177, 208, 213, 259 of globular cluster 
M3 from B06 and 8 variables 1, 6, 11, 13, 16, 18, 20, 24 of 
globular cluster M68 from W94 were used. We call this data 
set as B06+W94. In Fig. 10, we show the V-I inter-relations 
of the Fourier phase parameters (j>2i , </>3i , (j>4i for RRab and 
RRc stars respectively. Points lying above a threshold sigma 
(cthres) from the fit are rejected using multi-pass fitting al- 
gorithm in IDL . The threshold sigma level is chosen based 
on a reasonably high degree of correlation between the pa- 
rameters in the I and V band so that the number of data 
points for the fit retained are significant. The mean differ- 
ence of 0.61% in absolute magnitude, 0.23% in log(L/L0), 
4.03% in [Fe/H], 0.10% in effective temperature and 1.18% 
in mass will result when DF99 inter-relations are used. 

The Fourier inter-relations of the phase parameters 
from the I band to the V band have the following form: 



a + f3(f>n 



(7) 
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Figure 10. Observational inter-relations between V and I band Fourier phase parameters. The left and the right panels show the inter- 
relations for the RRab and RRc stars respectively. The solid line is the linear fit to the data. The inter-relations are based on the cosine 
series fit. Dots represent the data points from B06 and encircled dots from W94 . 



Table 3. Coefficients of the Fourier inter-relations for RRab stars 
from DF99 and B06+W94 data. 



Table 4. Coefficients of the Fourier inter-relations for RRc stars 
from DF99 and B06+W94 data. 



Relation 



4>2l' 



■ 031 

■ 041 



DF99 



B06+W94 





a 








P 


V 


0.693 


0.697 


2.097 ±0.152 


0.425 ±0.035 


V 


-0.039 


0.788 


0.436 ±0.075 


0.568 ±0.030 


V 


-0.611 


0.810 


0.868 ±0.132 


0.568 ±0.033 



Relat 


on 


DF99 


B06+W94 


021 J 
031 J 
041 ' 


-021 1 ' 
-031 V 
-041^ 


OL 

0.144 
-0.249 
-0.305 


P 
0.930 
0.995 
0.980 


a 

0.779 ± 0.302 
-0.083 ±0.081 
0.339 ±0.153 




0.803 ±0.062 
0.925 ±0.023 
0.672 ±0.069 



In Tables 3 and 4, we list the inter-relations obtained 
from the B06+W94 dataset for RRab and RRc stars. Also 
given are the theoretical inter-relations of DF99. In Fig. 10, 
we plot the I and V band Fourier parameters. In Fig. 11, we 
plot the theoretical inter-relations between the I and V band 
of DF99 and observational inter-relations obtained from the 
B06+W94 dataset. 

To compute V band amplitude from I band amplitude 



we again resort to the existing work in the literature. In 
Fig. 12 we plot V band amplitudes as functions of I band 
amplitudes for RRL stars in 4 globular clusters, M68 from 
Walker (1994), IC 4499 from Walker & Nemec (1996), NGC 
1851 from Walker (1998) and M3 from Benko et al. (2006). 
Solid line denotes the linear relation as determined by a least 
square fit given by 



A v = 0.071 (±0.019) + 1.500 (±0.040) Ai 



(8) 
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Figure 11. Theoretical inter-relations of DF99 between V and I band Fourier phase parameters in comparison with the observational 
relations calculated from the combined data set of B06 and W94 (B06+W94). The inter-relations are based on the cosine series fit. 



5.2 Mean magnitudes in V band 

Calibration to the standard V band magnitude is necessary 
for calculation of the distance scale. The mean magnitude 
was calculated for each object from its light curve. To trans- 
form the I band mean magnitude to the V band we used RRL 
light curves in the V and I band from the OGLE database 
itself. Typical errors in the V band magnitudes in OGLE 
are 0.02 - 0.07 mag. For RRab stars a linear least-square fit 
resulted in the following conversion relation : 

A (V) = -0.691 (±0.029) + 1.068 (±0.002) A (I). (9) 

For the RRc stars the result of the linear least square fit is 
as follows: 

A (V) = 3.733 (±0.536) ± 0.831 (±0.028) A (I). (10) 

The fits of Eqns. (9) & (10) are shown in Fig. 14. Light curves 
having mean magnitude variation ^ 2a were discarded and 
V band magnitude thus determined. On the other hand, the 
absolute magnitude calculation of RRab stars required Ai 
while RRc stars required A4 in V band. We used V and I 
band Data from Walker (1994, 1998) for the A 4 V band cali- 



bration and a linear regression analysis yielded the following 
relation 

Ai(V) = -0.007 (±0.001) + 1.686 (±0.040) Ai(I), (11) 
for RRab stars and 

A 4 (V) = 0.001 (±0.0007) ± 1.056 (±0.115) A 4 (I), (12) 

for RRc stars. Following Saha & Hoessel (1990) and Sakai 
et al. (1999), the intensity weighted mean magnitudes and 
phase-weighted mean magnitudes are given by 

n 

m int = -2.51og^il0- a4m ', (13) 

i=l 
n 

m ph = -2.5 log ^ 0.5(^+1 -0i-i)lO- a4mi , (14) 

i=l 

where n is the total number of observations, mi and (pi are 
the magnitude and phase of the ith observation respectively, 
in order of increasing phase. The intensity-weighted mean 
magnitude and phase-weighted mean magnitude in the V 
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Figure 12. V band amplitudes as functions of I band amplitudes 
for RRL stars in 4 globular clusters. Solid line denotes the linear 
relation as determined by a least squares fit. The diamonds denote 
the observational data taken from Walker 1994 (M68), Walker 
& Nemec 1996 (IC 4499) and Walker 1998 (NGC 1851). The 
triangles denote the amplitude relations determined by us from 
the published light curves from Benko et al. (2006) (M3). The 
quadratic correlation coefficient is 0.923. 
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Figure 13. Period-amplitude diagram for all RRab (open circles) 
and RRc stars (filled circles). 

band have been determined using a linear least square fit 
between the V and I band data resulting in the following 
conversion relations 

m int (V) = 0.862 (±0.360) + 0.984 (±0.019) m int (I), (15) 
m ph (V) = 0.471 (±0.320) + 1.009 (±0.021) m ph (I), (16) 



for RRab stars, and 

m int (V) = 3.093 (±0.867) + 0.863 (±0.045) m int (I), (17) 

m ph (V) = 5.900 (±0.893) + 0.733 (±0.045) m ph (I), (18) 
for RRc stars. 



6 PHYSICAL PARAMETER EXTRACTION 
6.1 RRab stars 

Using a set of 81 RRab stars, JK96 derived an empirical 
relation between the Fourier decomposition parameter </>3i, 
period (P) and metallicity [Fe/H] given by 

[Fe/H] = -5.038 - 5.394P + 1.345 03i. (19) 

The Fourier parameters in the above relation were obtained 
for a sine Fourier series fit. Care has to be taken that the 
RRab stars to which this equation is applied follow the light 
curve systematics defined by calibrating the sample of JK96. 
This can be ensured if the light curve in question satisfies 
a certain compatibility criterion. The various parameters of 
the calibrating data set obey well-defined correlations (see 
Table 6 of JK96). Deviation of a light curve from the cal- 
ibrating characteristics is quantified via the deviation pa- 
rameter Df defined as 

F obs - F calc 
Df = , 

OF 

where F b s is the observed value of a given Fourier param- 
eter, of is the corresponding deviations of various correla- 
tions. If the maximum value of Df of each parameter is less 
than 3, then each light curve satisfies the compatibility con- 
dition. This maximum value of D m represents a quality test 
on the regular behaviour of the shape of the light curve. But 
when the test is applied to the RRab stars of OGLE SMC 
RR lyrae stars, some of the stars did not pass this criterion 
even though the shapes of the light curve are 'normal look- 
ing'. One of the possible reasons is that the OGLE SMC 
RRL light curves are constructed from the observation in 
the I band but not in the V band on which the deviation 
parameter is defined. We have noticed that for the deviation 
parameter up to a value of D m = 5, the shapes of the RRab 
light curves are normal. This value of D m is used as a cutoff 
limit for the selection of clean samples of RRab stars for 
analysis. The deviation parameter of individual RRab light 
curves of the OGLE data has been calculated and is given 
in the last column of Table 1. Cacciari et al. (2005) have 
also adopted this slightly relaxed criterion to improve the 
statistics after they verified that this does not lead to any 
significant difference in the resulting physical parameters. 
For the RRab analysis, we apply various empirical relations 
from the literature to all our target stars only when D m ^ 5. 
After implementing this compatibility test we have found 
that 355 RRab stars out of 478 are now available for further 
analysis. This means that our sample of RRab stars now 
contains 355 RRab stars that have 'normal looking' light 
curve shapes. It is also seen that the maximum contribution 
to the deviation parameter comes from the deviation in </>3i . 

The intrinsic colors as derived from the Fourier parame- 
ters of the RRab stars can be used to estimate the tempera- 
ture of these stars. These color indices as defined by Jurcsik 
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Figure 14. Mean magnitude V band calibration for RRab stars (open circles) and RRc stars (filled circles). Top panels show fits 
represented by Eqns. (9) & (10). The lower panels show mean magnitude calibration using both V and I band data available for 465 
RRab and 55 RRc stars in SZ02. 



(1998) are the differences of the magnitude-averaged abso- 
lute brightness. The intrinsic color relations is of the form 



(B — V) = 0.308 + 0.163P- 0.187 Ai. 



(20) 



The above relation can be used to calculate the effective 
temperature 

logT cff (B - V) = 3.930 - 0.322 (B - V) + 0.007 [Fe/H]. (21) 

The relation between the absolute magnitudes and the 
Fourier parameters is given by the following relations 



M v = 1.221 - 1.396P - 0.477 Ai +0.103 031. 



(22) 



Eqn. (20) is from Jurcsik (1998) and Eqn. (21) is taken 
from Kovacs & Walker (2001). Eqn. (22), taken from Kovacs 
(1998) is based on the relation between the intensity- 
averaged My and the V band Fourier parameters. On the 
other hand on a different absolute scale, Jurcsik (1998) de- 
rived a relation for the mass of the RRab stars as follows 



log (M/M ) = 20.884 - 1.754 log P + 1.477 
log (L/L Q ) - 6.272 log T cff + 0.0367 [Fe/H]. 



(23) 



The Fourier coefficients in the above equations are 
based on a sine series fit. But the coefficients in the present 



analysis (Table 1) are calculated for a cosine series fit. There- 
fore, to put the coefficient 03i on to the appropriate system, 
we add 3.145 to 4>3i m Table 1. In Table 5 we list various 
physical parameters computed from the above empirical re- 
lations for 335 RRab stars with regular light curves with 
D m ^ 5 . Mean error in the Fourier phase parameter </>3i for 
335 RRab stars is ~ 0.19. This should result in an uncer- 
tainty of ~ 0.25 dex in [Fe/H], ~ 0.02 mag in the absolute 
magnitude, ~ 12 K in temperature, ~ 0.01 Mq in mass. 



6.2 RRc stars 

For the RRc stars, we have only retained stars for which 
the error in tfisi is less than 0.5. It is generally not desirable 
to have error larger than 0.2 (Kaluzny et al. 1998). Since 
the RRc stars in the SMC are fainter and have lower ampli- 
tudes, good data for such sources are simply not available at 
present. Imposing the criterion of <t$ 31 < 0.3 leaves us with 
only 3 RRc stars to be used for further analysis. To increase 
the statistics without significant loss of accuracy, we have 
chosen cr^, 31 < 0.5 which allows us to analyse 17 RRc stars. 
Through Hydrodynamical pulsation models, SC93 for- 
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mulated the following relationships for Luminosity, mass and 
helium abundance (Y) of RRc stars : 

log (L/L ) = 1.04 log P - 0.058 3i + 2.41, (24) 
log (M/M Q ) = 0.52 log P - 0.11 3 i + 0.39, (25) 
log (T off ) = 3.7746 - 0.1452 log P + 0.0056 3i . (26) 

log Y = -20.26 + 4.935 log (T cff ) - 0.2638 

log (M/M Q ) + 0.3318 log (L/Lq). (27) 

The rms error to the fit in Eqn. (24) is 0.025. 

Another alternative to Eqn. (24) is the expression for 
the intensity averaged absolute magnitude proposed by 
Kovacs (1998, hereafter K98) based on observations of RRc 
variables in eight different stellar systems: 

M v = 1.261 - 0.961 P - 0.044 021 - 4.447 A 4 , (28) 

with an rms of 0.042. Kovacs derived this formula from the 
light curves of 106 RRc stars, mainly from Sculptor and M68 
and calibrated the zero point from the Baade-Wesselink lu- 
minosity scale of Clementini et al. (1995). The phase differ- 
ence 021 in Eqn. (28) is based on a sine series fit. Therefore, 
1.571 is subtracted from 02i values in Table 2, which are 
based on a cosine series fit. 

Following Lazaro et al. (2006), we have calculated the 
luminosities of SMC RRc stars from Eqn. (24) as well as 
from Eqn. (28). The results are plotted in Fig. 15. Lumi- 
nosities calculated using SC93 relation (Eqn. (24)) are larger 
than those estimated by the K98 relation (Eqn. (28)) and the 
difference may be attributed to different scaling (zero points) 
of the two relations. The zero point of the luminosity scale 
of RRc stars is still a matter of controversy (Cacciari et al. 
2005, Nemec 2004, Lazaro et al. 2006). Some problems with 
the Simon's calibration for the RRc stars were also noted by 
Catelan (2004) who concluded that SC93's relations for L 
and M cannot both be simultaneously valid. Eqns. of SC93 
(Eqns. (24) & (25)) yield 

log P = -2.877 + 1.305 log (L/L ) -0.688 log (M/M Q ), (29) 

The above relation lacks a temperature term which 
should be present according to the period-density relation 
(Catelan 2004). This is true and remains an unresolved the- 
oretical issue. 

The empirical relation for metallicity in terms of the 
Fourier decomposition parameter 3 i was derived by Mor- 
gan et al. (2007) by considering a large number of stars in 
different globular clusters which have the V band Fourier co- 
efficients available in the literature and the sources of metal- 
licities available for these cluster members. Using the metal- 
licty scales of Zinn & West (1984) and Carretta & Gratton 
(1997), they obtained the following two empirical relations 
for RRc stars respectively 

[Fe/H] = 52.466 P 2 - 30.075 P + 0.131 3 i 2 + 

0.982 03i - 4.198 P03i + 2.424, (30) 

which has a sample standard deviation of 0.145 dex, and 

[Fe/H] = 0.034 3 i 2 + 0.196 3 i - 8.507 P + 0.367, (31) 



which has a sample standard deviation of 0.142 dex. We 
have used Eqn. (31) for the calculation of [Fe/H]. 

The temperatures determined using Eqn. (26) are found 
to be overestimated for some of the RRc stars. Cacciari et 
al. (2005) also used Eqn. (26) and found a similar trend 
for temperature estimates of RRc stars of globular cluster 
M3 when compared with temperatures derived by Sekiguchi 
& Fukugita (2000) from the B-V Colors. Temperature dif- 
ferences as large as 500 K have been found for some of the 
RRc stars. Therefore, the temperature of RRc stars as deter- 
mined from Fourier coefficients need to be re-examined. The 
inaccuracy of the empirical relations connecting the effective 
temperature, 03i term and the period is also reflected in the 
calculation of the radius of RRc stars. We show ahead that 
the radii of RRc stars determined from the Fourier decom- 
position technique differ largely from those determined from 
the theoretical period-radius-metallicity relation of Marconi 
et al. (2005). On the other hand, for RRab stars there is a 
well-matched similarity between the radii determined from 
the Fourier decomposition technique and from the more fun- 
damental relations used by Marconi et al. (2005). Therefore, 
the empirical relations derived for the RRab stars seem to 
be accurate, whereas the physical parameter estimation us- 
ing the empirical relations derived for the RRc stars may be 
unreliable. In Table 6, we list the various physical param- 
eters computed from the above empirical relations for the 
RRc stars with <t$ 31 < 0.5. 

Mean error in the Fourier phase parameter 03i for 17 
RRc stars is ~ 0.396. This will result in an uncertainty of 
~ 0.14 dex in [Fe/H], ~ 0.02 in log (L/L ), - 16 K in tem- 
perature and ~ 0.18 M in mass. 

6.3 HR Diagram 

An H-R diagram of the SMC RRLs stars is plotted in Fig. 15 
along with the theoretical blue and red edges of the insta- 
bility strip. The RRc stars are represented by solid circles 
whereas the RRab stars are represented by open circles. Also 
shown are the blue and red edges of the instability strip of 
Fundamental (RRab) and first-overtone (RRc) variables de- 
rived from the theoretical convective pulsation models of 
Bono, Caputo & Marconi (1995). It can be easily seen that 
all the RRab and RRc stars lie well within their instabil- 
ity strips. Although the RRc stars lie inside the instability 
strip, the empirical relations overestimate the temperatures 
of these stars. 

6.4 Gravity for RRab and RRc stars 

From the global physical parameters like mass, luminosity 
and temperature, the gravity can be calculated from the 
following equation (cf. Cacciari et al. 2005) 

log g = -10.607 + log (M/M ) - log (L/L ) + 4 log T eff (32) 



6.5 Radii of RRLs 

The radii of the RRL can be obtained once we know the 
temperature and the luminosity of the RRLs stars using the 
Stefan-Boltzmann law. Marconi et al. (2005) have recently 
given a new theoretical period-radius-metallicity relation for 
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Figure 15. H-R diagram of 335 RRab and 17 RRc stars and 
the theoretical pulsational blue and red edges of the instability 
strip from Bono et al. (1995). Open circles represent RRab stars, 
solid circles represent RRc stars with log(L/LQ) calculated from 
Eqn. (28) and filled upper triangles RRc stars with log(L/Lo) 
calculated from Eqn. (24). The green solid and dashed lines show 
the boundaries of the instability strip for RRc and RRab stars 
respectively for mass 0.75Mq and the black solid and dashed lines 
show the boundaries of the instability strip for RRc and RRab 
stars respectively for mass O.65M0. 



the RR Lyraes based on detailed and homogeneous set of 
nonlinear models with a wide range of stellar masses and 
chemical compositions. We show how the radii obtained by 
the Fourier decomposition method compare to the radii ob- 
tained from the period-radius-luminosity (PRZ) relation of 
Marconi et al. (2005). For the sake of completeness, their 
PRZ relations are given as follows: 

log R = 0.774(±0.009) + 0.580(±0.007) log P 

-0.035(±0.001) log Z, (33) 

for RRab stars with a = 0.008. 

log R = 0.883(±0.004) + 0.621(±0.004) log P 

-0.0302(±0.001) log Z, (34) 

for RRc stars with a — 0.004. log Z can be calculated from 
the relation 



log Z = [Fe/H] - 1.70 + log (0.638 f+ 0.362), 



(35) 



where f is an a-enhancement factor with respect to iron 
(Salaris, Chieffi & Straniero 1993). We take f = 1. In Fig. 16, 
we compare our determinations of the radii using the Fourier 
decomposition technique (FD) with those obtained by Mar- 
coni et al. (2005) empirical relations (Mar). In case of the 
radius determinations of RRc stars by the FD method, the 
SC93's equations for log L and log T are used. It can be 
seen that the radii obtained by the two different methods 
are quite consistent for RRab stars. The solid line in Fig. 16 
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Figure 16. Radii obtained from the empirical relation of Mar- 
coni et al. (2005) are plotted versus the radii determined from 
the Fourier decomposition (FD) method. For RRab (open cir- 
cles) stars in the upper panel the linear regression analysis 
yields : log (R/R©) (Mar) = -0.112 (±0.005) + 1.206(±0.007) log 
(R/R©) (FD) and for the RRc stars (filled circles) in the lower 
panel the linear regression analysis yields: log (R/R©) (Mar) = 
0.235(±0.064)±0.651(±0.090)log (R/R©) (FD). 



represents the best least square fit. The correlation coeffi- 
cient between the two radii determinations for the RRab 
stars are 0.994, whereas for the RRc stars (lower panel of 
Fig. 16), the correlation coefficient is 0.880. This suggests 
that the radii determinations by the two independent meth- 
ods are strongly correlated for the RRab stars, whereas this 
is not the case for the RRc stars. 



7 DISTANCE TO THE SMC 

We use mean magnitude Ao(V), intensity- weighted mean 
magnitude and phase weighted mean magnitude separately 
to derive an independent distance modulus to the SMC. The 
distance modulus for each star is calculated from their val- 
ues of mean magnitudes and My derived from the Fourier 
parameters. Assuming the reddening estimates of the eleven 
SMC fields of SZ02 to be accurate, we take average of the 
reddening values of these fields as the true estimates of E(B- 
V). Following SZ02 we take the interstellar extinction as 
Ay-=3.24E(B-V). This value of A v has been used to esti- 
mate the reddening free distance modulus of SMC. Using 335 
RRab and 17 RRc stars we find the mean distance modulii 
of SMC to be 18.89±0.01 mag, 18.87± 0.03 mag and 18.87 
±0.02 mag from the mean magnitude, intensity-weighted 
mean magnitude and phase-weighted mean magnitude re- 
spectively for all the 352 RRL stars. The uncertainty is the 
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standard deviation of the mean from the individual stars. 
The mean distance modulus so determined are consistent 
with the earlier values of 19.05± 0.017 (Kovacs et al. 2000), 
18.89±0.4 (Harries et al. 2003) and 18.91±0.1 (Hilditch et 
al. 2005). 



8 SUMMARY & CONCLUSIONS 

In this paper we have derived the physical parameters of 
352 RR Lyrae stars (335 RRab and 17 RRc) of the SMC 
from OGLE-II I band database using the Fourier decompo- 
sition of their light curves. The stars were selected based 
on the quality of their light curves. I band Fourier coef- 
ficients have been converted to V band using the inter- 
relations obtained from the observational data of B06 and 
W94. Using the Fourier decomposition method, we find the 
mean physical parameters: [Fe/H] = -1.56 ±0.25, M = 0.55 
±0.01 M s , T eff = 6404 ± 12 K, logL = 1.60 ± 0.01L o 
and M v = 0.78 ± 0.02 for 335 RRab variables and [Fe/H] 
= -1.90 ± 0.13, M = 0.82 ±O.18M , T cff = 7177 ± 16 K, 
log L = 1.62 ± 0.02 L Q and M v = 0.76 ± 0.05 for 17 RRc 
stars. Mean distance modulus to the SMC was calculated 
from the 352 light curves by using mean magnitude, inten- 
sity weighted mean magnitude and phase weighted mean 
magnitude. The values of the distance modulus are found 
to be in good agreement with independent studies. Loca- 
tions of the RRL stars in the H-R diagram clearly show 
that the estimates of the parameters determined from the 
Fourier decomposition method are consistent with the theo- 
retical blue and red edges of the instability strip calculated 
by Bono et al. (1995). The calculations of the radii of the 
RRLs by using the Fourier decomposition technique are in 
agreement with the theoretical period-radius-metallicity re- 
lations of Marconi et al. (2005) obtained from a completely 
different approach of non-linear convective modelling. 
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Table 5. Physical parameters extracted from Fourier coefficients for 335 RRab variables. Errors represent the uncertainties in the Fourier 
parameters. 
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Complete table is available in the electronic form. 

Table 6. Physical parameters extracted from Fourier coefficients for 17 RRc variables. Errors represent the uncertainties in the Fourier 
parameters. 



OGLE ID M V (K98) log(L/L )(SC93) [Fe/H] T eS M/M logg R/R (FD) R/R (Mar) 





(1) 






(2 


) 




(3) 




(4) 






(5) 




(6) 






(7) 




(8) 




(9) 


005451 


.72-723850. 


4 





.85± 


0.03 


1 


66± 


0.02 


-1 


,10± 


0. 


17 


7457±17 





59± 





19 


3 


.00± 


0.33 


4. 


07± 


0.03 


4. 


19± 


0.02 


005115 


.64-724739. 


2 





,78± 


0.04 


1 


67± 


0.02 


-1 


21± 





16 


7425±17 





62± 





10 


3 


.00± 


0.31 


4. 


19± 


0.04 


4 


25± 


0.02 


010245 


.99-721132. 


7 





.83± 


0.04 


1 


83± 


0.03 


-2 


02± 


0. 


11 


7165±18 


1 


13± 





07 


3 


,04± 


0.06 


5 


35± 


0.05 


4 


63± 


0.02 


004631 


.02-724658. 


8 





,78± 


0.04 


1 


,75± 


0.03 


-1 


,67± 


0. 


17 


7266±20 





77± 





.18 


2 


,97± 


0.24 


4 


79± 


0.05 


4 


65± 


0.02 


004902 


.13-724513. 


6 





.82± 


0.03 


1 


76± 


0.02 


-1 


73± 





13 


7247±15 





80± 





.18 


2 


.97± 


0.23 


4. 


87± 


0.04 


4 


69± 


0.02 


010453 


.59-723756. 








,78± 


0.03 


1 


83± 


0.01 


-2 


,11± 


0. 


OS 


7115±10 





93± 





17 


2 


.94± 


0.18 


5 


44± 


0.04 


5 


09± 


0.01 


010029 


.57-725454. 


2 





,77± 


0.04 


1 


,81± 


0.02 


-2 


03± 


0. 


00 


7136±10 





85± 





10 


2 


.93± 


0.22 


5 


32± 


0.05 


5 


10± 


0.01 


010029 


.57-725454. 


2 





,77± 


0.04 


1 


,81± 


0.02 


-2 


03± 





00 


7136±10 





85± 





10 


2 


.93± 


0.22 


5 


32± 


0.05 


5 


10± 


0.01 


004327 


.25-724542. 


4 





,76± 


0.03 


1 


,75± 


0.02 


-1 


62± 





13 


7238±13 





64± 





22 


2 


.89± 


0.34 


4. 


80± 


0.04 


4. 


96± 


0.02 


010231 


17-722134. 








69± 


0.04 


1 


83± 


0.03 


-2 


,13± 





15 


7107±18 





91± 





.18 


2 


.93± 


0.19 


5 


46± 


0.06 


5 


16± 


0.02 


005556 


.74-732133. 


6 





,77± 


0.05 


1 


,77± 


0.03 


-1 


80± 





17 


7193±18 





70± 





21 


2 


.89± 


0.31 


5 


00± 


0.06 


5 


07± 


0.03 


003934 


.35-730433. 








,76± 


0.03 


1 


79± 


0.02 


-1 


92± 





13 


7161±13 





76± 





21 


2 


.90± 


0.28 


5 


16± 


0.04 


5 


13± 


0.02 


010647 


.29-723053. 


7 





,73± 


0.04 


1 


,81± 


0.02 


-2 


03± 





13 


7124±14 





.77± 





21 


2 


.88± 


0.27 


5 


31± 


0.04 


5 


27± 


0.02 


005155 


.37-724909. 


5 





73± 


0.04 


1 


,77± 


0.02 


-1 


78± 





17 


7183±16 





65± 





23 


2 


.86± 


0.35 


5 


00± 


0.05 


5 


19± 


0.03 


010316 


.37-724816. 


2 





66± 


0.12 


1 


90± 


0.03 


-2 


56± 





13 


6968±18 


1 


14± 





12 


2 


.92± 


0.11 


6 


20± 


0.15 


5 


54± 


0.02 


010316 


.37-724816. 


2 





66± 


0.12 


1 


90± 


0.03 


-2 


56± 





13 


6968±18 


1 


14± 





12 


2 


.92± 


0.11 


6 


20± 


0.15 


5 


54± 


0.02 


005527 


.97-724136. 


8 





,76± 


0.04 


1 


81± 


0.03 


-2 


03± 





Hi 


7118±17 





75± 





22 


2 


.87± 


0.29 


5 


32± 


0.05 


5 


35± 


0.03 



Kovacs, G., Kupi, G. 2007, A&A, 462, 1007 

Lazaro, C, Arellano Ferro, A., Arevalo, M. J. et al. 2006, 

MNRAS, 372, 69 

Marconi, M., Nordgren, T., Bono, G. et al. 2005, ApJ, 623, 
L133 

Morgan, S. M., Simet, M., Bargenquast, 1998, Acta As- 
tron., 48, 341 (MSB98) 

Morgan, S. M., Wahl, J. N., Wieckhorst, R. M., 2006, Mm- 
SAI, 77, 178 

Morgan, S. M., Wahl, J. N., Wieckhorst, R. M., 2007, MN- 
RAS, 374, 1421 

Nemec, J., M., 2004, AJ, 127, 2185 

Pefia, J. H., Arellano, A., Sareyan, J. P., Pena, R., Alvarez, 



M. 2007, Comm. in Asteroseismology, 150, 385 
Petersen, J. O., 1986, A&A, 170, 59 

Press, W., Teukolsky, S., Vetterling, W., Flannery, B. 1992, 
Numerical Recipes in FORTRAN (2nd ed.; Cambridge: 
Cambridge Univ. Press) 

Pritzl, B. J., Armandroff, T. E., Jacoby, G. H., Da Costa, 

G. S., 2005, AJ, 129, 2232 

Saha, A., Hoessel, J. G., 1990, AJ, 99, 97 

Sakai, S., Ferrarese, L., Kennicutt, R. C. et al., 1999, ApJ, 

523, 540 

Salaris, M., Chieffi, A., Straniero, O., 1993, ApJ, 350, 645 
Sandage, A., 2004, AJ, 128, 858 

Schaltenbrand, R., Tammann, G., 1971, A&AS, 4, 265 



Physical parameters of the SMC RR Lyrae stars 



Sekiguchi, M., Fukugita, M., 2000, AJ, 120, 1072 

Simon, N. R., Clement, C. M., 1993, ApJ, 410, 526 (SC93) 

Simon, N., Lee, A. 1981, ApJ, 248, 291 

Smith, H. A., Silbermann, N. A., Baird, S., R., Graham, J. 

A., 1992, AJ, 104, 1430 

Soszyriski, I., Udalski, A., Szymahski, M. et al., 2002, Acta 

Astron., 52, 369 (SZ02) 

Walker A. R., 1994, AJ, 108, 555 (W94) 

Walker A. R., Nemec, J. M., 1996, AJ, 112, 2026 

Walker A. R., 1998, AJ, 116, 220 

Walker, A. R., Mack, P., 1988, AJ, 96, 872 

Zinn, R., West, M. J., 1984, ApJS, 55, 45 



